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Using a cleanly tagged data sample of charged current events, it is demonstrated that the rate 
at which such events are mis-identified as v e 's is accurately simulated in the MiniBooNE — » v E 
analysis. Such mis-identification, which could arise from muon internal bremsstrahlung, is decisively 
ruled out as a source of the low energy electron-like events reported in the MiniBooNE search for 
Vn —* Ve oscillations. This refutes the conclusions of a recent paper which postulates that hard 
bremsstrahlung could form a substantial background to the MiniBooNE v e sample. 



The MiniBooNE Collaboration has reported the re- 
sults of a search for — > v e oscillations at Am 2 ~ l 
eV 2 . In this search, no significant excess of events was ob- 
served above background for reconstructed neutrino en- 
ergies above 475 MeV, but 96 ± YJ{stat) ±20(sys)_ excess 
events were reported between 300 and 475 MeV [lj . The 
data are not consistent with two-neutrino oscillations and 
the source of the excess is, at present, undetermined. 

A recent paper @ suggests that the source of the low 
energy excess is muon internal bremsstrahlung associated 
with the i/fj, charged current quasi-elastic (CCQE) inter- 
action, i/fj, n fi~ p + 7, as depicted in Fig[T] The is- 
sue of bremsstrahlung of muon neutrino CCQE events 
was raised in an earlier paper Q and brought to Mini- 
BooNE 's attention in the course of the oscillation analysis 
prompting the study presented here. 

Reference asserts that v^n — > pT p+7 events can be 
mis-classified as v e CCQE signal events if the final state 



muon is below Cerenkov threshold, since the radiated 
photon may imitate a final state electron. Even when the 
muon escapes direct detection, however, its presence may 
be revealed by the Michel electron from the muon decay. 
Making use of this Michel electron tag, the present pa- 
per demonstrates that the misidentification of events 
caused by muon internal bremsstrahlung cannot be the 




FIG. 1: Feynman diagram for muon internal bremsstrahlung. 
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FIG. 2: Reconstructed neutrino energy distribution of one 
subevent events passing u e selection cuts. The data are shown 
as points with statistical errors and the Monte Carlo predic- 
tion is shown as a histogram with systematic errors. This 
figure is a reproduction of Fig. 2 in 



source of the observed low energy excess and is not a 
significant background to the MiniBooNE — ► v e oscil- 
lation search. Rather than relying on Monte Carlo simu- 
lation of this process, data are used to directly constrain 
the contribution of the muon bremsstrahlung diagram. 
The study was conducted prior to the unblinding of the 
data for the oscillation analysis and the results were in- 
corporated into the estimated backgrounds at that time. 
Because the study showed that the background to the os- 
cillation analysis from muon internal bremsstrahlung was 
extremely small, it was deemed unnecessary to add such 
internal radiative effects to the simulation. The analysis 
makes use of an event sample in which the presence of 
a muon is tagged strictly by the presence of the Michel 
electron from the decay of the muon: 



V +P, M 



The MiniBooNE detector and trigger [l[ are particu- 
larly well suited for such identification. The MiniBooNE 
trigger creates a 19.2 [is time window surrounding the 1.6 
/xs beam spill. When events are reconstructed, periods of 
time in which light is produced in the tank are identified 
as "subevents". Subevents are separated by looking for 
gaps between PMT hit times larger than 10 ns, and are 
typically ^100 ns in length. 

In 82% of the cases where a muon is contained in the 
detector, a second subevent from the Michel electron is 
produced. The 18% without a second subevent splits 
into 8% that result from fi~ capture in oil (this rate has 
been separately measured Q), 2% where simulation pre- 
dicts the Michel electron creates too few PMT hits to 
be clearly seen (<10 hit PMTs), and 8% where the muon 
decays sufficiently quickly that the decay cannot be time- 
resolved from the initial interaction. 

To study the rate at which charged current interac- 
tions arc mis- identified as f e 's, events with two subevents 
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FIG. 3: Reconstructed neutrino energy distribution of two 
subevent events passing v e selection cuts. The data are shown 
as points with statistical errors and the Monte Carlo pre- 
diction is shown as a histogram. No systematic uncertainty 
has been evaluated for this sample. Normalization of both 
data and Monte Carlo is to the rate of CCQE events with 
no observed Michel. If muon internal bremsstrahlung were 
the source of the low energy excess in the oscillation sample 
(Fig. HJ , then the same excess would also be observed in this 
figure, but is not. 



are first selected. The two subevents must be separated 
by at least 1500 ns to avoid any instrumental effects 
causing the second subevent. No cuts are placed on the 
first subevent, but the second subevent must have fewer 
than 6 veto hits to ensure containment and reject cosmic 
muons and fewer than 200 main tank hits to ensure an 
energy consistent with a Michel electron. Having created, 
in this way, a sample of events tagged as being from a 
Up, charged current interaction, the second subevent and 
its hits are discarded. The full set of v e selection cuts, 
identical to those used in the oscillation result, are then 
applied to this artificial one subevent sample. The re- 
sulting sample is a direct, in situ measurement of the 
bulk of the charged current contribution to the back- 
ground in the v e sample, including that from muon inter- 
nal bremsstrahlung. The only i/ M charged current back- 
ground component this measurement misses is that due 
to muon decay so rapid that the Michel electron cannot 
be time separated from the parent muon. This back- 
ground is constrained and checked in other ways. 

The v e selection cuts [l| include precuts to isolate a 
clean neutrino event sample: the event must have just 
one subevent within the 1.6 /is beam window and have 
fewer than 6 veto hits to remove incoming cosmic rays 
and exiting muons from neutrino events in the detector. 
Each event must also have more than 200 tank hits to re- 
ject Michel electrons from stopped cosmic rays which can 
enter the tank prior to the trigger window and therefore 
avoid the veto. 

After the precuts, the same "track-based" algorithm 
used in the oscillation analysis [l[ reconstructs the vertex 
position, angle, energy, and time of the event, assuming 
the light comes from an extended, straight line source. 
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The vertex and projected track endpoint must lie within 
the fiducial volume of the detector and the visible energy 
in the tank must be E vis > 140 McV. 

Events are identified as v e CCQE based on likelihoods 
which are built from phototube charge and time prob- 
ability distribution functions. For each event, using a 
single-track hypothesis, the likelihood is calculated that 
it is an electron (L e ) or a muon (L^). The event is then 
reconstructed under a two track hypothesis, where the 
invariant mass is forced to be 135 MeV, and a tt° likeli- 
hood {L v ) is formed. Finally, the event is reconstructed 
by finding the best two track fit allowing the invariant 
mass to float. This yields a best- fit mass (M 77 ). Visible 
energy-dependent cuts on \og(L e J L^), \og(L e / L^) and 
M 77 are then applied to isolate a v e CCQE signal sam- 
ple. The neutrino energy is then determined using the 
reconstructed lepton energy and direction and assuming 
the interaction is v e CCQE. 

Fig. [2] shows the spectrum of the originally-published 
v e candidate sample subject to this selection The ex- 
cess of events between 300 and 475 MeV is clearly visible. 
Fig. [3] shows the results of the same cuts applied to the 
tagged data sample. In this case, the normalization 
has been adjusted to correct for the Monte Carlo deter- 
mined rates of one subevent and two subevent i/ M charged 
current interactions and for the requirement that the two 
subevents be separated by more than 1500 ns. Fig. [3] is 
therefore a direct prediction and measurement of the 
charged current contribution to Fig. [2] The Monte Carlo 
prediction for the rate at which these processes pass the 
v e selection agrees well with the data. Furthermore, the 
abscissa of these two figures shows that charged cur- 
rent processes account for only a tiny fraction of the back- 
ground in Fig. [2] Muon internal brcmsstrahlung (or any 
Up, charged current process that can lead to muon decay 



at rest) therefore cannot be the source of the low energy 
data excess. 

Muon brcmsstrahlung events should be efficiently re- 
jected by the v e selection. Muon tracks are considerably 
longer than electron tracks for the same visible energy in 
the MiniBooNE detector. For events where the muon en- 
ergy lies above Cerenkov threshold, the additional track 
length leads to a significantly different charge and hit 
structure, particularly in the center of the Cerenkov ring. 
One would therefore expect the presence of the muon to 
pull the charge and time likelihoods of the events away 
from an "electron-like" hypothesis. 

In this paper, charged current events are identified 
through the presence of a Michel electron and then sub- 
jected to the MiniBooNE oscillation analysis cuts. The 
rate of misidentification of these events as ^ e 's is accu- 
rately modeled by the Monte Carlo and is not the source 
of the low energy excess in MiniBooNE. 

It should be noted that this study says nothing about 
the total rate of muon internal brcmsstrahlung in Mini- 
BooNE, just the rate at which this process occurs in the 
v e sample. It may be possible to relax or adjust the v e 
selection cuts to make a measurement of muon internal 
brcmsstrahlung at some time in the future. 
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